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From sand piles to electrons in metals, one of 
the greatest challenges in modern physics is to un- 
derstand the behavior of an ensemble of strongly 
interacting particles. A class of quantum many- 
body systems such as neutron matter and cold 
Fermi gases share the same universal thermo- 
dynamic properties when interactions reach the 
maximum effective value allowed by quantum me- 
chanics, the so-called unitary limit [l|, Qj. It 
is then possible to simulate some astrophysical 
phenomena inside the highly controlled environ- 
ment of an atomic physics laboratory. Previous 
work on the thermodynamics of a two-component 
Fermi gas led to thermodynamic quantities aver- 
aged over the trap [H-dj], making it difficult to 
compare with many-body theories developed for 
uniform gases. Here we develop a general method 
that provides for the first time the equation of 
state of a uniform gas, as well as a detailed com- 
parison with existing theories [6-15]. The preci- 
sion of our equation of state leads to new phys- 
ical insights on the unitary gas. For the unpo- 
larized gas, we show that the low-temperature 
thermodynamics of the strongly interacting nor- 
mal phase is well described by Fermi liquid the- 
ory and we localize the superfluid transition. For 
a spin-polarized system [16-18], our equation of 
state at zero temperature has a 2% accuracy and 
it extends the work of (l9l . I20I ] on the phase dia- 
gram to a new regime of precision. We show in 
particular that, despite strong interactions, the 
normal phase behaves as a mixture of two ideal 
gases: a Fermi gas of bare majority atoms and a 
non-interacting gas of dressed qua si-particles, the 
fermionic polarons 



01 dressed quasi 

[10, flSE-Si]. 



In this letter we study the thermodynamics of a mix- 
ture of the two lowest spin states (i = 1,2) of 6 Li pre- 
pared at a magnetic field B = 834 G (see Methods), 
where the dimensionless number l/kpa characterizing 
the s-wave interaction is equal to zero, the unitary limit. 
kp is the Fermi momentum and a the scattering length. 
Understanding the universal thermodynamics at unitar- 
ity is a challenge for many-body theories because of the 
strong interactions between particles. Despite this com- 
plexity at the microscopic scale, all the macroscopic prop- 
erties of an homogeneous system are encapsulated within 
a single equation of state P(/zi, H2, T) that relates the 




FIG. 1: Schematic representation of the universal function 
h(n,C)- It fully describes the thermodynamics of the unitary 

function of chemical potential imbalance r\ = 
and of the inverse of the fugacity £ = exp(— /ii/fcsT). In 
this paper we measure the function h over the black lines 
(1, Q and (77, 0) which correspond to the balanced unitary gas 
at finite temperature and to the spin-imbalanced gas at zero 
temperature, respectively. 



pressure P of the gas to the chemical potentials Hi of the 
species i and to the temperature T. In the unitary limit, 
this relationship can be expressed as 



(1) 

is the pres- 



P(fii,(i 2 ,T) = Pi(ni,T)h ( rj = — ,( = exp 

V Mi 

where Pi( M i,T) - -k B TX^(T)f 5/2 (-C" 1 ) 
sure of a single component non-interacting Fermi gas and 
/ 5 / 2 (z) = X)^i z "/ n5 ^ 2 - hiViQ i s a universal function 
which contains all the thermodynamic information of the 
unitary gas (Fig. 1). In cold atomic systems, the inhomo- 
geneity due to the trapping potential makes the measure- 
ment of h(rj, C) challenging. However, this inhomogeneity 
of the trap can be turned into an advantage as shown in 

We directly probe the local pressure of the trapped 
gas using in situ images, following the recent proposal 
[23j . In the local density approximation, the gas is locally 
homogeneous with local chemical potentials: 



Mi (r) = Ml ° - V(r) 



(2) 



is the chemical potential at the bottom of the trap 
for species i). Then a simple formula relates the pressure 
P to the doubly- integrated density profiles [23] ]: 
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FIG. 2: Schematic representation of our atomic sample. The 
6 Li atomic cloud is imaged in the direction y; the column 
density is then integrated along the direction x to give n(z). 
The 7 Li atoms are imaged after a time of flight along the z 
direction. 



where rii(z) = J rii(x, y, z)dxdy, rii being the atomic 
density. uj t (resp. ui z ) is the transverse (resp. axial) 
angular frequency of a cylindrically symmetric trap (see 
Fig. 2) and \n z = m(0,0,z) is the local chemical po- 
tential along the z axis. From a single image, we thus 
measure the equation of state (JTJ) along the parametric 
line (77, C) = (p2z/mx,exp(-Vu/kBT)), see below. 

The interest of this method is straightforward. First, 
one directly measures the equation of state (EOS) of the 
uniform gas. Second, each pixel row Zi gives a point 
h(r](zi)X( z i)) whose signal to noise ratio is essentially 
given by the one of rii(z) + 712(2); typically one experi- 
mental run leads to ~ 100 points with a signal to noise 
between 3 and 10. With about 40 images one gets ~ 4000 
points h(rj, C), which after averaging provides a low-noise 
EOS of standard deviation a = 2%. In the following 
we illustrate the efficiency of our method on two impor- 
tant sectors of the parameter space (t),£) in Fig. 1: the 
balanced gas at finite temperature (1,0 an d the zero- 
temperature imbalanced gas (77, 0). 

We first measure the equation of state of the unpo- 
larized unitary gas at finite temperature, P{jjL\, n%,T) — 
P(/j,,T). The measurement of h(l,() through the local 
pressure ([3]) can be done provided one knows the temper- 
ature T of the cloud and its central chemical potential fi . 

In the balanced case, mo del- independent thermometry 
is notoriously difficult because of the strong interactions. 
Inspired by [24| , we overcome this issue by measuring the 
temperature of a 7 Li cloud in thermal equilibrium with 
the 6 Li mixture (see Methods). 

/x° is fitted on the hottest clouds so that the EOS agrees 
in the classical regime £ 3> 1 with the second-order virial 
expansion h(l, C) ^ 2(l+C~ 1 /'V / 2) [2j|. For colder clouds 
we proceed recursively. The EOS of an image recorded 
at temperature T has some overlap with the previously 
determined EOS from all images with T" > T. In this 
overlap region fjr is fitted to minimize the distance be- 
tween the two EOS's. This provides a new portion of the 
EOS at lower temperature. Using 40 images of clouds 
prepared at different temperatures, we thus reconstruct 
a low-noise EOS from the classical part down to the de- 
generate regime, as shown in Fig. 3a. 

We now comment the main features of the equation of 



state. At high temperature, the EOS can be expanded 
in powers of C 1 as a virial expansion (llj : 

ft(i,C)_££i ((-i) k+1 k-^ + b k )c- k 



2 ET=i(- 1 ) k+lk ~ 5/2 (~ k 

where b^ is the k th virial coefficient. Since we have 62 = 
l/y/2 in the measurement scheme described above, our 
data provides for the first time the experimental values of 
63 and 64. 63 = —0.35(2) is in excellent agreement with 
the recent calculation 63 = — 0.291 — 3 -5 / 2 = —0.355 from 
11] but not with b 3 = 1.05 from [ll]. 64 = 0.096(15) 
involves the 4-fermion problem at unitarity and could 
interestingly be computed along the lines of 

Let us now focus on the low-temperature regime of the 
normal phase (< 1. As shown in Fig. 3b, we observe 
a T 2 dependence of the pressure with temperature. This 
behavior is reminiscent of a Fermi liquid and indicates 
that pseudogap effects expected for strongly-interacting 
Fermi superfluids [26| do not show up at the thermody- 
namic level within our experimental precision. In analogy 
with 3 He or heavy-fermion metals, we fit our data with 
the EOS: 



P(M,T) = 2P 1 ( M ,0) & 
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PifaO) = l/157r 2 (2m//i 2 ) 3 / 2 /i 5 / 2 being the pressure of 
a single-component Fermi gas at zero temperature, to* 
is the quasi-particle mass and £~ 1 is the compressibility 
of the normal gas extrapolated to zero temperature, and 
normalized to that of an ideal gas of same density. We 
deduce two new parameters m*/m = 1.13(3) and £ n = 
0.51(2). Despite the strong interactions to* is close to to, 
unlike the weakly interacting 3 He liquid for which 2.7 < 
TO*/m < 5.8, depending on pressure. Our £ n value is in 
agreement with the variational Fixed-Node Monte-Carlo 
calculations £„ = 0.54 in 27 1, £„ = 0.56 in [lj| and with 
the Quantum Monte-Carlo calculation £„ = 0.52 in [28| . 
This yields the Landau parameters Fg = ^„to*/to — 1 = 
-0.42 and F( = 3(to*/to - 1) = 0.39. 

In the lowest temperature points (Fig. 3c) we ob- 
serve a sudden deviation of the data from the fit (0| at 
(kBT/n) c = 0.32(3) (see supplementary materials). We 
interpret this behavior as the transition from the normal 
phase to the superfluid phase. This critical ratio has been 
extensively calculated in the recent years. Our value is 
in close agreement with the diagrammatic Monte- Carlo 
calculation (^r/fi) c = 0.32(2) of and with the Quan- 
tum Monte-Carlo calculation {k B T/n) c = 0.35(3) of [HI 
but differs from the self-consistent approach in Q giv- 
ing (hsT / 'fj,) c = 0.41, from the renormalization group 
prediction 0.24 in [2^], and from several other less pre- 
cise theories. From eq. (j4]) we deduce the total density 
n = ni + 712 — dP((j,i — /^,T)/9/u and the Fermi energy 
Ep = ksTp = fr 2 /2m(3n 2 n) 2 / 3 at the transition point. 
We obtain {ji/E F ) c = 0.49(2) and (T/Tf) c =0.157(15), 
in very good agreement with [6j. Our measurement is 
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the first direct determination of (n/Ep) c and (T/Tp) c 
in the homogeneous gas. It agrees with the extrapolated 
value of the MIT measurement [13] ■ 

Below T c , advanced theories 0, S] predict that 
P(fi,T)/2P l ((j l ,0) is nearly constant (Fig. 3b). Therefore 
at T = T c , P/2Pi ~ £7 3/2 ^ 3.7, and is consistent with 
our data. Here £ s = 0.42(1) is the fundamental parame- 
ter characterizing the EOS of the balanced superfluid at 
zero temperature, a quantity extensively measured and 
computed in the recent years 0. 

Our data is compared at all temperatures with the cal- 
culations from (Fig. 3a). The agreement with 0] is 
very good for a large range of temperatures. Concern- 
ing [6[, the deviation with our data is about one error 
bar of the Monte-Carlo method below £ = 0.2 and the 
deviation increases with temperature (Fig. 3a). Fur- 
thermore, we show in the supplementary material that 
/i(l,C)/2 must be greater than 1, an inequality violated 
by the two hottest Monte-Carlo points of Q. 

From our homogeneous EOS we can deduce the equa- 
tion of state of the harmonically trapped unitary gas 
by integrating h(l,Q over the trap (see supplementary 
material). In particular, we find a critical tempera- 
ture for the trapped gas (T/Tp) c = 0.19(2), where 
Tp = h^uiluizN) 1 / 3 . This value agrees very well with 
the recent measurement of 30], and with less precise 
measurements @, [3ll, 

Let us now explore a second line in the universal dia- 
gram h(r], £) (Fig. 1) by considering the case of the T = 
spin-imbalanced mixture (12 7^ Ml, i-S- f] 7^ 1- Previous 
work [IB - [l8j ] has shown that phase separation occurs in 
a trap. Below a critical population imbalance a fully- 
paired supernuid occupies the center of the trap. It is 
surrounded by a normal mixed phase and an outer rim 
consisting of an ideal gas of the majority component. In 
two out of the three previous experiments including ours 
[l6L [Tij , the local density approximation has been care- 
fully checked. We are therefore entitled to use ([3]) to 
analyze our data. 

As in the previous case, the relationship between the 
pressure and the EOS requires the knowledge of the 
chemical potentials /i° and ^§ a t the center of the trap. 

/ii is determined using the outer shell of the ma- 
jority spin component (i = 1). The pressure profile 
P (fJ'iz, A*2z, 0) corresponds to the Fermi-Dirac distribu- 
tion and is fitted with the Thomas-Fermi formula Pi 



oil - z 2 /P 2 ) 5/2 , providing /x? = 



\R{. Using Pi 



for the calculation of h = PjP\ cancels many system- 
atic effects on the absolute value of the pressure. More- 
over, fitting the outer shell using a finite-temperature 
Thomas- Fermi profile (l9| . we measure a temperature 
k B T = 0.03(3)^?. 

^2 is fitted by comparison in the superfluid region with 
the superfluid equation of state at zero temperature [2lj : 

h( V ,0) = (l + v f/y(2i s ) 3 /\ (5) 

Our measured equation of state h(rj, 0) is displayed in 
Fig. 4. By construction our data agrees for 77 > 0.1 with 
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FIG. 3: (color online) Equation of state of a spin-balanced 
unitary Fermi gas. (a) Finite-temperature equation of state 
h(l, C) (black dots). The error bars represented at £ = 0.14 
and £ = 2.3 indicate the 6% accuracy in £ and h of our EOS. 
The red curves are the successive virial expansions up to 4 th 
order. The blue triangles are from the green stars from FS, 
the purple diamonds from [8[ , and the blue solid line from . 
The grey region indicates the superfluid phase, (b) Equation 
of state P(/i,T)/2Pi(/i,0) as a function of (k B T/n) 2 , fitted 
by the Fermi liquid equation of state ((4]). The red dashed 
line is the non- interacting Fermi gas (NIFG). The horizontal 
dot-dashed (resp. dotted) line indicates the zero-temperature 
pressure of the superfluid phase oc S^ 3 ^ 2 (resp. normal phase 
oc £n )• (c) Expanded view of (b) near T c . The sudden de- 
viation of the data from the fit occurs at (fc_BT/ j u) c = 0.32(3) 
that we interpret as the superfluid transition. The black 
dashed line indicates the mean value of the data points below 



eq.©. In Fig. 4 the slope of h(r),0) displays an obvious 
discontinuity for 77 = r\ c = 0.065(20). This is a signature 
of a first-order quantum phase transition to the partially 
polarized normal phase. The error bar is dominated by 
the uncertainty on £ s . This value is slightly higher than 
the prediction r\ c — 0.02 given by the fixed-node Monte- 
Carlo [l0( and than the value r\ c = 0.03(2) measured in 

From the relations rii — dPj d^t we deduce from h(i~i, 0) 
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FIG. 4: (color online) Equation of state of the zero- 
temperature spin-imbalanced unitary gas h{n,Q) (black dots). 
Error bars are equal to one standard error. The red solid line 
is the superfiuid equation of state, the blue dashed line is the 
ideal Fermi liquid equation with A = —0.615, m* = 1.20m 
and the black dotted line is the Monte Carlo calculation from 
[Tot . Inset: Local density ratio n%/ni as a function of r\. The 
red solid line 712/111 = 1 corresponds to the fully paired su- 
perfiuid and blue dashed line to the model @. 



the density ratio n^jnx (see inset in Fig. 4). This ratio is 
discontinuous at the phase transition, from a maximum 
value in the normal phase (ri2/n,i) c — 0.5(1) to n% 
11 in the superfiuid phase. Our value is close to the 
zero-temperature calculation 0.44 [lOj and agrees with 
the coldest MIT samples [l|J [2(|. It confirms that the 
temperature is much smaller than the tricritical point 
temperature T = 0.07Tp [lj| where the discontinuity 
vanishes, justifying our T — assumption made above. 

For r] < rj c our data displays a good agreement with 
a simple polaron model, based on the pioneering work 
in [ldj . A polaron is a quasi-particle describing a sin- 
gle minority atom immersed in the majority Fermi sea 
flol 18, 21, [2^. It is characterized by a renormalized 
chemical potential fj,2 — A fi\ and an effective mass m* 
10]. Following this picture, we write the pressure as 



the sum of the Fermi pressures of ideal gases of majority 
atoms and of polarons: 



P 



1 / 2m 



3/2 



mP + (5) 3/2 ^-^) 5/2 ], 

(6) 
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which can be written as 



fm*\ 3/2 
%0) = 1+Mj (n-A) 
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(7) 



A and m* have recently been calculated exactly [lj, LL5J : 
A — —0.615, m*/m = 1.20(2) and with these values in- 
serted in (O the agreement with our data is perfect. Note 
that our data lies slightly above the variational fixed- 
node Monte Carlo calculation [lOj. We therefore con- 



b 3 64 {k B T/n) c (u/E F ) c (T/2» c 

-0.35(2) 0.096(15) 0.32(3) 0.49(2) 0.157(15) 

£ n m* jm T} c (n 2 /ni) c m*/m 

0.51(2) 1.13(3) 0.065(20) 0.5(1) 1.20(2) 

TABLE I: Table of quantities measured in this work. 



elude that interactions between polarons are not visible 
at this level of precision. 

Alternatively, we can fit our data with m*/m as a free 
parameter in ©. We obtain m*/m = 1.20(2). The 
uncertainty combines the standard error of the fit and 
the uncertainty on £ s . This value agrees with our pre- 
vious measurement m*/m = 1.17(10) [l8| (with a 5- 
fold improvement in precision) , with the theoretical value 
m*/m = 1.20(2) in 14, TEJ and with the variational cal- 
culation Hp It differs from the values 1.09(2) in [33| . 
1.04(3) in [10[ , and from the experimental value 1.06 in 
0. 

We arrive at a simple physical picture of the T = 
spin-polarized gas: the fully paired superfiuid is de- 
scribed by an ideal gas EOS renormalized by a single 
coefficient £ s ; the normal phase is nothing but two ideal 
gases, one of bare majority particles and one of polaronic 
quasi-particles. 

In conclusion, we have introduced a powerful method 
for the measurement of the equation of state of the uni- 
tary and homogeneous Fermi gas, that enables direct 
comparison with theoretical models and provides a set 
of new parameters shown in Tab[l] The method can 
readily be extended to any multi-component cold atom 
gas in three dimensions that fulfills the local density ap- 
proximation (see supplementary discussion). We have 
shown that the normal phase of the unitary Fermi gas is a 
strongly correlated system whose thermodynamic prop- 
erties are well described by Fermi liquid theory, unlike 
high-r c cuprates. 

Note added in proof: Since this paper was accepted Re- 
publication, we have become aware of the measurement of 
a similar equation of state for the balanced unitary Fermi 
gas at finite temperature by different methods [34 1. 
METHODS SUMMARY 

Our experimental setup is presented in 18j . We load into 



an optical dipole trap and evaporate a mixture of 6 Li in 
the 1 1/2, ±1/2} states and of 7 Li in the |1,1) state at 
834 G. The cloud typically contains N 6 = 5 to 10 x 10 4 
6 Li atoms in each spin state and Nj = 3 to 20 x 10 3 7 Li 
atoms at a temperature from T = 150 nK to 1.3 /iK. The 
6 Li trap frequencies are u z /2ir = 37 Hz, uj r /2ir varying 
from 830 Hz to 2.20 kHz, and the trap depth is 25 fiK for 
our hottest samples, with T ~ 2Tp. 6 Li atoms are im- 
aged in situ using absorption imaging, while 7 Li atoms 
are imaged after time of flight, providing the tempera- 
ture in the same experimental run (Fig. 4). Since the 
scattering length describing the interaction between 7 Li 
and 6 Li atoms, = 2 nm, is much smaller than hp 1 , 
the 7 Li thermometer has no influence on the 6 Li density 
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profiles. The 7 Li- 6 Li collision rate, Tq? = 10 s _1 , is large 
enough to ensure thermal equilibrium between the two 
species. 
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METHODS 

Construction of the equation of state by succes- 
sive patches. A typical image at high temperature 
provides about 100 pixels corresponding to £ values 
varying from 2 at the trap center to 6 at the edges, with 
a signal-to-noise from 3 to 10. 7 such images are fitted 
in the wings using the second-order virial expansion and 
averaged to obtain a low- noise EOS up to £ = 2. Then 
images of clouds where the evaporation has been pushed 
to a slightly lower temperature are recorded. They show 
about 75% overlap in ( with the previous EOS. After 
minimization of the distance between a new image and 
the previously determined EOS in the overlap region, 
we obtain the value of /x° for a single image with 3% 
statistical uncertainty. This process is repeated for 6 
successive trap depths. When averaging one image with 
typically 10 previous images, we obtain a new EOS 
with an error on C of about 0.03/^ ~ 1%. The EOS 
experiences a random walk error on the 40 images of 
0.01 x \/~40 ~ 5% for the coldest data. An independent 
check of the maximum error is provided by the good 
agreement with the superfluid equation of state for 
temperatures lower than T c 0, H| . 

Evaluation of the systematic uncertainties. For 

the measurement of /i(l,£), the combined uncertainties 
on the radial frequency of the trap, trap anharmonicity, 
magnification of our imaging system, and atom counting 
affect the pressure measurement given in Q at ~ 20% 
level. However, two measurements, one at relatively high 
temperature and one at very low temperature, enable us 
to show that the overall error does not exceed 6%. In 
the temperature range £ > 0.5, the agreement between 
the experimental value b$ = —0.35(2) and the theoretical 
value 63 = —0.355 of the third virial coefficient indicates 
that the global systematic error is smaller than 6%. Sec- 
ond, at very low temperature, theory 0, Q predicts that 
the variation of Pj2P\ as a function of fcsT//i in the su- 
perfluid phase remains smaller than 5%. Our value of 
P/2Pi = 3.75 below the critical point is within 5% of 
the T = prediction £7 3/2 = 3.7(2). This confirms that 
systematic errors for our coldest samples are also smaller 
than 6%. 

For the determination of the critical transition 
to superfluidity we fit the low-temperature data 
P(/x, T)/2Pi(/Lt, 0) with a variable horizontal line for T < 
T c and with the Fermi-liquid equation (j4|) for T > T c . 
The result of the fit is the dashed black line in Fig. 3c, 
which intersects equation (j4]) at (fcsT/^) c = 0.315(8). 
This statistical error is negligible compared to the er- 
ror induced by the 6% systematic uncertainty discussed 
above, justifying our very simplified fit procedure. In- 
deed a 6% error on the pressure induces a 10% error on 



/i for images recorded in the vicinity of the critical tem- 
perature, leading to (fcsT//^) c = 0.32(3). 

For the measurement of h(r], 0), the fit of the fully po- 
larized wings of the cloud serves as a pressure calibration 
for the rest of the cloud, cancelling many systematic ef- 
fects. 

In order to estimate temperature effects in the polar- 
ized gas, let us first remark that in the superfluid phase 
corrections scale as T 4 for the bosonic excitations and 
are exponentially suppressed by the gap for the fcrmionic 
ones [7|. So in our temperature range fc^T = 0.03 fi® their 
contributions will be very small. On the other hand, in 
the partially polarized normal phase, we expect a typical 
Fermi liquid T 2 scaling. In order to obtain an estimate 
of the error on the EOS, we develop the following sim- 
ple model. In equation which describes a mixture of 
zero-temperature ideal gases, we replace the Fermi pres- 
sures by the finite-temperature pressures of ideal gases 
(see equation ((TJ): 

fm*\ 3/2 

P(H 1> H 2) T)=P 1 (fx 1 ,T)+ (-£) P 1 (fi 2 -A f i 1 ,T), 

and run the analysis described in the main text. At T = 
0.05/i?, the correction on h is less than 1%, half of our 
current error bar. 

Limit of 7 Li Thermometry. As the scattering length 
between the 7 Li atoms, 077 = —3 nm is negative, the 
7 Li cloud becomes unstable when a BEC forms. This oc- 
curs at T ~ 150 nK with typically 3500 atoms. Precise 
thermometry with lower atom numbers becomes difficult. 
For the measurement of the zero-temperature equation of 
state of the imbalanced gas, we do not use 7 Li thermom- 
etry but rather the fit of the wings of the majority spin 
component. 
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SUPPLEMENTARY DISCUSSION 
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FIG. 5: (color online) Equations of state of the trapped uni- 
tary gas. (a) Comparison between our EOS E/NEf as a 
function of S/Nks (black dots) and the EOS measured in 
[fj (open red squares). The red solid line is the second-order 
virial equation of state, (b) Comparison between our EOS 
E/E° as a function of (T/T F )° (black dots) and the EOS 
measured on 40 K in || (open red squares). The grey regions 
correspond to the superfluid phase. 



Equation of State of the Trapped Unitary Gas 

In this work, we have measured the equation of state 
of the homogeneous unitary gas. We can deduce from 
our data the EOS of the trapped balanced unitary gas, 
which has been measured in fl, HI . 

Using the local density approximation, the total atom 
number N = J n dr 3 is expressed as a function of the 
temperature T and the chemical potential /j, at the cen- 
ter, involving the function h(l,Q: 



N 



dlog^C/Co) 



i-k \ huj ) J Co dC 

(8) 

where Co = exp(— pP/ksT) and uj = (lo^ujz) 1 / 3 . We 
use for the calculation a discretized version of (j8]) taken 
solely on our experimental values of h, i.e. without us- 
ing any interpolating or fitting function. Similar ex- 
pressions are used to calculate the Fermi temperature 
E F = k B T F = huj^N) 1 / 3 , the total entropy S and en- 
ergy E of the cloud. The equation of state E/NE F as a 
function of S/Nk B , displayed in Supplementary FigJSK, 
is in very good agreement with 

The normal-supernuid phase transition for the trapped 
gas occurs when at the trap center Co = Cc = 
exp(-(fc B T/^)- 1 ), with {k B T/n) c = 0.32(3), as mea- 
sured on the homogeneous EOS /i(l,£). At this point 
we get (T/T F ) C = 0.19(2), (S/Nk B ) c = 1.5(1) and 
(E/NE F ) C = 0.67(5). 

In order to make the comparison with Q, we also 
express the equation of state E/E° as a function of 
(T/T F )°, where the superscript refers to the quan- 
tities evaluated on a non-interacting Fermi gas having 



the same entropy (Supplementary FigfSJj). The good 
agreement with the measurement in [3j, performed on 
40 K clouds, illustrates the universality of the unitary gas. 



Physical interpretation of the pressure in the normal 
mixed phase 

We have shown that the pressure in the normal mixed 
phase can be described as the sum of the Fermi pressures 
of ideal gases of majority atoms and of polarons: 



P 




Am) 



5/2 



(9) 

Here, we evaluate the corresponding canonical EOS re- 
lating the energy E to the densities ni,ri2, in order to 
compare with the Fixed Node Monte Carlo prediction 
[ioj . Since at unitarity we have E = 3PV/2 we just 
express the chemical potentials in terms of densities by 
using the thermoc 
yield respectively: 



using the thermodynamical identities rii = 9 Mi P, which 



n-2 



ni 



1 

6^ 



2m* 



(fx 2 - A^) 



.3/2 



(10) 

(11) 



The last term in equation (fTTj) clearly indicates the in- 
creased majority density due to the presence of the mi- 
nority component. Expressing the pressure as a function 
of rii in (j9]) yields the energy: 



E = E, 



FP 



K 771* 



where E FF is the energy of the fully polarized gas and 
x = n^/nx. Expanding E to order x 2 finally leads to an 
expression similar to that obtained in 



E{x) = E FP 1 



-Ax 



r 5 / 3 



Bx 2 



with B = 5A 2 /9 = 0.2. Our value of B is close to the 
calculated value B ~ 0.14 from floj [ . 



Trap Anharmonicity 

First, in the axial direction z, the confinement is pro- 
duced magnetically and the corresponding anharmonic- 
ity is negligible. In the radial direction, we develop the 
gaussian potential to fourth order around p = 0: 



V r (p) = V 1-exp^- 



1 2 2 i 4 

-mu) r p + ep , 
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where muj = 2Vo/a 2 and e — — Vb/2cr 4 . In the balanced 
case, we have 



n(z) = J d 2 pn ffjP - i 



2 2 2 2 4 

niw'z — -mLJ r p — ep 



Introducing n = dP/dp and defining u = moj 2 p 2 /2 + ep 4 
we obtain, to lowest order, 



recovering the true scattering length, one has trivially 
(V)q < 0, hence f2 < Oq. Using the thermodynamic 
identity ft — —PV, and recalling that flo = —2P{V and 
h = P/Pi, we finally get the inequality 

Mi»0>2. 



2tt 



The error on the measurement of h is then 



V 



muj 2 n(z) ^ _ k B T 



2wP 1 ( f x z ,T) 



/s/aK'-^Ml.C) , w 



v J c / 5/ 2(-C" 1 ) C 



(12) 

We evaluate the integral in (fT2"j) using the experimental 
values of h(l,Q- hi our shallowest trap, the worst case 
anharmonicity effect is 5%. 



Extension to a Multi-Component System 

We extend the equation (2) to a mixture of species 
i, of mass m,-, trapped in a harmonic trap of transverse 
frequencies w r i, following the calculations in 23]. Us- 
ing Gibbs-Duhem relation at a constant temperature T, 
dP = J2i n i dp l , then 



E 



2tt 



E 



2tt 



■dxdy 



<9P 



/E d ^ 



<9P 

dpi 



An exact inequality on the equation of state of an 
attractive Fermi gas 

Writing the hamiltonian as H = Hq + U, where Hq is 
the single-particle part of the hamiltonian and U is the 
inter-particle interaction, one has the general inequality 
il < flo + (V)o, where fio is the grand potential associated 
with Hq and (-)o is the thermal average related to Hq [35j | . 
Taking for U a short range square potential of depth Uq < 



where we have used local density approximation [pi (r) = 
p® — V(r)) to convert the integral over space to an inte- 
gral on the chemical potentials. The integral is straight- 
forward and yields to 

P(p lz ,T) = — ^2 m^riiiz). 



